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Abstract. Realistic 3-dimensional (3D), radiative hydro- 
dynamical surface convection simulations of the metal- 
poor halo stars HD 140283 and HD 84937 have been per- 
formed. Due to the dominance of adiabatic cooling over 
radiative heating very low atmospheric temperatures are 
encountered. The lack of spectral lines in these metal- 
poor stars thus causes much steeper temperature gradi- 
ents than in classical ID hydrostatic model atmospheres 
where the temperature of the optically thin layers is deter- 
mined by radiative equilibrium. The modified atmospheric 
structures cause changes in the emergent stellar spectra. 
In particular, the primordial Li abundances may have been 
overestimated by 0.2-0.35 dex with ID model atmospheres. 
However, we caution that our result assumes local thermo- 
dynamic equilibrium (LTE), while the steep temperature 
gradients may be prone to e.g. over-ionization. 
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1. Introduction 

Stellar chemical compositions are of great astrophysical 
importance as they carry information on stellar evolution 
and mixing processes, galaxy formation and evolution, and 
Big Bang nucleosynthesis and the baryon density of the 
Universe. In order to decipher the observed stellar spectra 
in terms of abundances a proper understanding of the line 
formation process is required. Often-used assumptions and 
approximations in abundance analyses involve ID, plane- 
parallel, hydrostatic and flux-constant (radiative and con- 
vective equilibrium) model atmospheres in LTE. Free pa- 
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rameters such as the mixing length parameters, micro- 
and macroturbulence enter the model constructions and 
the spectral syntheses, and cause uncertainties in the de- 
rived abundances. 

An improved theoretical foundation for the inter- 
pretation of stellar spectra is desirable, to complement 
the impressive recent observational advances in terms of 
signal-to-noise, spectral resolution and limiting magni- 
tudes. Modern supercomputers now allow self-consistent 
3D, radiation- hydrodynamics simulations of the surface 
convection of late-type stars. These simulations invoke no 
free parameters, yet succeed in reproducing key diagnos- 
tics such as granulation topology and statistics, helioseis- 
mic properties, and spectral line shapes, shifts and asym- 
metries for the best available test-bench: the Sun (Stein 
& Nordlund 1998; Rosenthal ct al. 1999; Asplund et al. 
1999b). Since convection may be expected to have a promi- 
nent influence on the atmospheres of metal-poor halo stars 
we have performed 3D simulations for a number of such 
stars. In the present Letter, we briefly present results for 
two of the stars and compare emergent spectra from the 
3D models with those of conventional ID models. A full 
account of these and other metal-poor convection simula- 
tions will be presented elsewhere. 



2. Hydrodynamical surface convection simulations 

Realistic 3D, time-dependent, surface convection sim- 
ulations of the metal-poor halo stars HD 140283 and 
HD 84937 have been performed using a compressible 
radiative-hydrodynamical code, that has previously been 
successfully applied to studies of the solar granulation 
(Stein & Nordlund 1998). The equations of mass, momen- 
tum and energy conservation coupled to the 3D equation 
of radiative transfer are solved on an Eulerian mesh with 
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Table 1. Adopted stellar parameters 



Star 


2tff 


logs 


[Fe/H] 


ta 

Sturb 




[K] 


[m s" 2 ] 




[kms" 1 ] 


HD 140283 


5690 


1.67 


-2.50 


1.0 


HD 84937 


6330 


2.04 


-2.25 


1.0 



a Not necessary for the 3D spectral line synthesis 



100 x 100 x 82 zones, covering 35 x 35 x 12 Mm (HD 140283) 
and 21x21x8 Mm (HD 84937), respectively. 

In order to accurately describe the photospheric lay- 
ers, special care must be taken to include appropriate in- 
put physics. A state-of-the-art equation-of-state (Mihalas 
et al. 1988), which includes the effects of ionization, exci- 
tation and dissociation, has been used together with de- 
tailed continuous (Gustafsson et al. 1975 and subsequent 
updates) and line (Kurucz 1993) opacities. The 3D radia- 
tive transfer is solved for one vertical and four inclined 
rays under the approximations of local thermodynamical 
equilibrium (LTE) and grouping of the opacities into four 
bins (Nordlund 1982). The accuracy of the opacity bin- 
ning procedure is verified throughout the simulations by 
solving the full monochromatic radiative transfer (2748 
wavelength points) in the 1.5D approximation and found 
to always agree within 1% in emergent flux. It is notewor- 
thy that the convection simulations contain no adjustable 
parameters besides those specifying the effective temper- 
ature T e ff (or, as used here, the entropy of the inflowing 
gas at the lower boundary) , the gravitational acceleration 
at the surface logg and the chemical composition [Fe/H]. 
The adopted stellar parameters, which have been obtained 
from the infrared flux method (IRFM, Alonso et al. 1996) 
for T c ff, Hipparcos parallaxes for logg and published val- 
ues for [Fe/H], are listed in Table |l|. The individual el- 
emental abundances are taken from Grevesse & Sauval 
(1998) scaled appropriately. Further details on the analo- 
gous solar simulations may be found in Stein & Nordlund 
(1998). 

Fig. [J shows the resulting temperature structure for 
the HD 140283 simulation, together with the correspond- 
ing solar simulation (Asplund et al. 1999b). The similarity 
between the two simulations in the deeper layers is partly 
fortuitous, and is due to an almost perfect cancellation of 
the effects of different log ,9 and [Fe/H]. More striking are 
the prominent differences in the temperature structure of 
the optically thin layers. While the solar simulation re- 
mains close to radiative equilibrium, much lower temper- 
atures are encountered in the HD 140283 simulation. This 
seems to be a generic feature as is also present in our 
simulations of HD 84937 and other low-metallicity stars. 
The large differences also remain when comparing on op- 
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Fig. 1. Average T(p) for HD 140283 (solid) compared with 
the Sun (dashed). Thin solid and dashed curves represent the 
corresponding extreme temperatures for a given density 
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Fig. 2. T(t) for 3D and ID models of HD 140283 and HD 84937 

tical depth scales as done in Fig. || where the spatial (on 
surfaces of equal optical depths) and temporal averages 
from the 3D simulations are presented together with the 
corresponding ID MARCS structures. 

The temperatures in the optically thin layers are pre- 
dominantly determined by two competing effects: adia- 
batic cooling due to the expansion of ascending material 
and radiative heating by spectral lines (radiative heat- 
ing occurs when spectral lines reabsorb some of the con- 
tinuum photons whose release cause the intense cooling 
near continuum optical depth unity). Thus, in dynamical 
atmospheres the temperature tends to become depressed 
below its radiative equilibrium value and lines act mainly 
as heating agents. This is opposite to the case when ra- 
diative equilibrium is enforced in theoretical ID model 
atmospheres. There, spectral lines cause surface cooling 
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(e.g. Mihalas 1978), hence leading to shallower tempera- 
ture gradients when the metallicity is diminished. In the 
convection simulations, the opposite happens: when de- 
creasing [Fe/H], fewer and weaker lines are available to 
absorb photons, allowing the adiabatic cooling to domi- 
nate more. As a result, balance between heating and cool- 
ing is achieved at lower surface temperatures, causing a 
steeper T(t). 

3. Impact on spectral line formation 

Using the convection simulations as model atmospheres 
the 3D radiative transfer was computed under the assump- 
tion of LTE for lines of special astrophysical interest, as 
listed in Table From the full simulations representa- 
tive sequences of one hour with snapshots every 2 min 
were selected for the line calculations; the emergent fluxes 
for the shorter intervals correspond to T e g = 5672 ± 3 K 
(HD 140283) and T cS = 6356 ± 8K (HD 84937), respec- 
tively, i.e. very close to the intended T c g . Prior to the spec- 
tral synthesis the simulations were interpolated to a finer 
vertical resolution (r max = 300). The radiative transfer 
equation was solved for a total of 29 directions (4 //-angles 
and 7 azimuthal angles plus the vertical) after which a 
disk- integration was carried out. Elemental abundances 
were derived from equivalent widths, to allow direct com- 
parison with ID models. It should be noted that no micro- 
or macroturbulence enter the calculations, as the convec- 
tive Doppler shifts and velocity gradients are sufficient 
to produce the observed line broadening (Asplund et al. 
1999b). To estimate the impact of 3D models, a differ- 
ential comparison with classical ID MARCS model atmo- 
spheres (Gustafsson et al. 1975 with subsequent updates) 
was carried out, adopting a microturbulence of 1.0 km s . 
The stellar parameters for the ID models are identical to 
those for the shorter convection simulation sequences. 

The reliability of the simulations may be assessed from 
a comparison between observed and predicted line asym- 
metries. For HD 140283 the theoretical bisectors agree 
very well with observations (Allende Prieto et al. 1999; 
Asplund et al. 1999a), e.g. in explaining the greater line 
asymmetries in this metal-poor star compared with solar 
metallicity comparison stars. The good concordance pro- 
vide additional confidence in the simulations. 

Due to differences in their sensitivity to the tempera- 
ture structure, lines are influenced differently by the dif- 
ferences between 3D convection simulations and ID model 
atmospheres, as is evident from Table ^. Lines of neutral 
minority species, low excitation transitions and stronger 
lines are formed in higher layers and thus feel the low tem- 
peratures there, making them stronger for a given abun- 
dance. We note for example that from the resonance line of 
Li I a 0.2-0.35 dex lower abundance is derived from the 3D 
models than from the ID models. On the other hand, high 
excitation transitions and lines of ionized elements tend to 
result in slightly larger abundances for the 3D models, as 



Table 2. The derived LTE abundances of 3D hydrodynamical 
simulations compared with ID hydrostatic model atmospheres 
(on the customary logarithmic scale where logeH = 12.00) 



Line HD 140283 HD 84937 



H 


3D 


lD a 


3D 


lD a 


Li 1 670.8 


1.78 


2.12 


2.08 


2.28 


Be ii 313.1 


-0.92 


-1.06 


-0.85 


-0.88 


B 1 249.8 


-0.41 


-0.21 






O 1 777.2 


7.20 


7.15 


7.32 


7.30 


K 1 769.9 


2.60 


2.74 






Ca 1 616.2 


3.62 


3.77 


4.22 


4.35 


Fel 6 


4.57(0.16) 


5.02 (0.17) 


5.10(0.12) 


5.37(0.11) 


Fen 6 


5.16 (0.10) 


5.08 (0.11) 


5.38 (0.02) 


5.35 (0.03) 


Ban 614.2 


-1.28 


-1.12 


-0.22 


-0.05 



a Adopting £turb = 1.0 km s 1 

b 10 Fel lines, 2 Fell lines. The dispersions are in parenthesis 

is the case for Be II, O I and Fe II. The larger effects for 
the HD 140283 simulation compared with HD 84937 may 
be attributed to the generally stronger lines in the former 
due to its lower T c g . Finally, we stress that the main lesson 
to draw from Table |^ is not the absolute 3D abundances 
but the relative abundances compared with differential ID 
model atmospheres. 

4. Discussion 

It is noteworthy that ionization balance of Fe is not ful- 
filled for the two 3D simulations, in contrast to the ID 
MARCS models, which suggests that either 1) the adopted 
stellar parameters are inappropriate, 2) the atmospheric 
temperatures have been underestimated, or 3) there are 
significant departures from LTE present in the form of a 
general over-ionization of Fe of 0.3-0.6 dex. The latter ef- 
fect is certainly not inconceivable, given the steep T(t) 
and the weak UV-blocking in these metal-poor environ- 
ments; a 3D non-LTE investigation is clearly of high pri- 
ority. One may speculate that if indeed Fe is over-ionized 
compared with LTE expectations, Li may also show a sim- 
ilar effect, in spite of the generally small non-LTE correc- 
tions in the ID case (Carlsson et al. 1994). The possibility 
of erroneous T e g s can most likely be refuted, since prelimi- 
nary calculations reveal only minor differences 50 K) in 
IRFM estimates between our 3D simulations and the ID 
MARCS models. The uncertainties in the stellar parame- 
ters (AT c ff = 80 K, Alogg = 0.15) only correspond to 
A[Fei] = 0.10 and A [Fen] = 0.08 and can thus probably 
not be blamed for the discrepancy either. 

One may suspect that the very low surface tempera- 
tures have indeed been underestimated, causing system- 
atic effects on the derived abundances. Support for such 
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a suspicion comes from an observed trend between abun- 
dance and excitation potential for the Fe I lines, that is 
present with both metal-poor simulations but slightly less 
pronounced and in opposite direction with the ID mod- 
els; unfortunately the sample of lines is still somewhat 
limited however. As discussed above, it is physically very 
plausible to have sub-radiative equilibrium temperatures 
in these metal-poor environments, but the magnitude of 
the effect may have been overestimated. One possible rea- 
son could be the neglect of Doppler shifts in the radiative 
energy transfer employed in the convection simulations. 
In the optically thin layers the temperature is largely con- 
trolled by only a few strong lines that are partially sat- 
urated. Velocity gradients will enable the lines to absorb 
unattenuated continuum photons instead, thereby caus- 
ing additional heating. Improved simulations addressing 
also this remaining refinement is planned. We note, how- 
ever, that our solar simulations already predict essentially 
perfect agreement of line shapes and asymmetries with 
observations (Asplund et al. 1999b), as well as agreement 
for a range of other diagnostics. 

Though the results presented here should be consid- 
ered preliminary, pending further investigations of depar- 
tures from LTE and the reality of the low surface temper- 
ature, it is still interesting to estimate possible effects on 
Big Bang nucleosynthesis and galactic chemical evolution. 

Our 3D results suggest that the primordial Li abun- 
dance may have been overestimated by rj 0.25 dex when 
using classical ID model atmospheres. Adopting log eu = 
2.20 ± 0.03 as the best estimate of the Li-plateau from a 
sample of metal-poor halo stars on the same T e ff-scale as 
used here (Bonafacio & Molaro 1997), we conclude that 
the primordial Li abundance may be as low as log era = 
1.95 ± 0.03, corresponding to Li/H=8.9 ± 0.6 ■ 1GT 11 . It 
should be noted that the la uncertainties only represent 
the (ID) statistical errors (Ryan et al. 1999) and no correc- 
tions for departures from LTE, post-Big Bang nucleosyn- 
thesis production of 6+7 Li or stellar Li-depletion have been 
applied; the uncertainties are undoubtedly dominated by 
systematic errors. Additional simulations of metal-poor 
stars are important, in order to estimate more carefully 
the statistical errors, and also to check whether the ex- 
treme thinness of the Li-plateau found when using ID 
models (Ryan et al. 1999) can be retained with 3D models. 
We note that our LTE Li estimate is only barely consis- 
tent with standard Big Bang nucleosynthesis predictions: 
if correct, our results imply that the baryon density in the 
Universe is pb ~ 1-7 • 10~ 31 gem -3 (77 w 2.5- 10 -10 ), which 
is in good agreement with 4 He (Y p = 0.238±0.002±0.005) 
but at variance with both of the contested high and low 
D measurements (e.g. Olive 1999). 

A discussion of the full impact of our results on galac- 
tic chemical evolution is clearly beyond the scope of the 
present paper. Instead, we very briefly point out some 
implications for the production of the light elements in 
the early Galaxy. Since Fe presumably suffers from over- 



ionization, [Fe] is best represented by [Feii3D] which is 
slightly higher than [Felio]- Therefore the slope [Be/H] 
vs [Fe/H] should remain essentially unchanged compared 
with ID, i.e. close to unity. Though OH lines have not been 
investigated here, we anticipate a decrease of the derived 
O abundance with the 3D simulations, and thus that the 
[Be/H] vs [0/H]-relation should be shallower than in ID. 
It would thus be more consistent with a primary origin 
for Be. Likewise we expect the [O/Fe] vs [Fe/H]-relation 
(based on OH-lines) to be flatter compared with the most 
recent findings (Israelian et al. 1998), but the claimed con- 
sistency between OH, 1 and [0 1] lines can no longer be 
expected, emphasizing the need for a non-LTE study of 
O based on these inhomogeneous atmospheres. Further- 
more, since Be 11 and B 1 react in opposite ways to the low 
atmospheric temperatures, the B/Be ratio may previously 
have been overestimated. However, departures from LTE 
driven by over-ionization are significant already in ID for 
B (Kiselman & Carlsson 1996) and it seems likely that 
in 3D the non-LTE effects will be aggravated, which may 
well be true also for other elements. 
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